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Disorder in the spectral function of a qubit ensemble 
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The impact of disorder in the vibrational bath of an ensemble of open quantum systems is explored, 
arising either from variation in the overall coupling strength or from uncertainty in the shape of 
the environment spectral function. Such disorder leads to an additional source of decoherence 
of subsystem ensembles, due to variation in the reorganization energy. Additionally, vibrational 
disorder induces a shift in the oscillation frequency of the coherence between ground and excited 
states. This shift is temperature dependent, and in the high temperature limit Suj oc T. This latter 
finding could be of particular relevance to biological systems exhibiting long lived coherences in 
highly disordered environments, where temperature-dependent quantum beats have been observed. 


I. INTRODUCTION 

The decay of quantum coherences is an active topic of 
research across a wide range of fields in physics, chem¬ 
istry and increasingly biochemistry n An important 
source of decoherence arises from interactions between 
the subsystem of interest and a vibrational heat bath 
The vibrational bath can be characterized by a 
continuous spectral density, J{'w), which quantifies the 
density of vibrational states with frequency w in the 
range 0 < w < oo, weighted by the strength of their cou¬ 
pling to the subsystem. While dissipative processes are 
often controlled by the finite-frequency structure of the 
spectral function, coherence losses due to pure dephasing 
are primarily controlled by the behaviour of this mec- 
tral function as it approaches zero frequency [1, [l,[a|. A 
high density of low frequency modes is expected to drive 
a rapid decay in coherences, while an absence of such 
modes may enable partial coherence to survive over long 
timescales. 

The pure dephasing properties of general, power-law 
spectral functions have been thoroughly investigated and 
have found important ^plications in quantum infor¬ 
mation and metrology [a, Simple pure dephasing 
models, such as the independent boson model, can be 
solved exactly Q. Recently, it was proposed that long- 
lasting quantum coherences observed in biological photo¬ 
synthetic systems could be stabilized by suppressed low 
frequency spectral densities, enabling long lived quantum 
coherences between excitonic states to survive at physio¬ 
logical temperatures. Kreisbeck and Kramer Q have in¬ 
vestigated this with advanced numerical methods for the 
Fenna-Matthews-Olson (FMO) complex, finding that the 
super-Ohmic behaviour of experimentally derived spec¬ 
tral functions for FMO greatly extend exciton-exciton co¬ 
herences without compromising efficient energy transfer 
0. A number of other theories, invoking highly struc¬ 
tured spectral densities typical of protein-bound chro- 
mophores, have also been presented 

Biological systems are known to exhibit a high degree 
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of static, structural and dynamic disorder [il[i3, yet 
typically their vibrational environments are modeled by 
well-defined spectral functions. The goal of this paper 
is to explore, within a simple model, the effect of in¬ 
corporating a random component to the mean spectral 
function. Though this problem is inspired by quantum 
biology, it has potential relevance across a much wider 
range of applications. 


II. A SIMPLE MODEL, THE TWO LEVEL 
SYSTEM 

Consider the independent boson model @,0, described 
by the Hamiltonian H = Hs -I- Hb + Hsb, and repre¬ 
sented schematically in figure [TJ 


Hs 

= |0 >< 0\+Ei\l >< 1|, 

(1) 

Hb 

k 

( 2 ) 

Hsb 

= ^{gkal^glak) |1 >< 1|. 

( 3 ) 


_ k 



FIG. 1: During direct excitation, the two energy levels are 
separated by an energy gap Ei. The vibrational states then 
respond; lowering the energy of the excited state to Ei — A, 
where A labels the reorganization energy. In the presence of 
vibrational disorder, a distribution of reorganization energies 
emerges, shown here with standard deviation cta. Realiza¬ 
tions of the system with smaller reorganization energies ex¬ 
perience less pure dephasing, leading to an apparent shift in 
oscillation frequency. 
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Hs represents a simple two level system. The up¬ 
per level is coupled to a set of vibrations, represented 
by bosonic oscillators. We can characterize the vibra¬ 
tional states by the continuous spectral density function 
= Y.k \9k\'^^{w - Wk) 0 . 

Imagine coherently exciting the system to the upper 
state at time t = 0. Since ground and excited states 
are decoupled, the population in the upper state cannot 
decay; this is represented by the diagonal component of 
the density matrix pii(t) = 1. At time zero the coherence 
between the ground and excited states will also be unity, 
represented by the off diagonal component poi(O) = 1- 
We are interested in the evolution of this coherence over 
time. 

Since the system was excited abruptly at t = 0, we 
are free to assume that the bath and the subsystem are 
initially unentangled; at t = 0 the bath lies in the thermal 
state /Tr[e~^^^]. Under such assumptions, the 

evolution of the system is known to obey 0 


Pl2{t) 

lit) 

fiw) 


Pi 2 ( 0 )e ' ""e 


-iEit^ 

dw J{w)f{w, t), 


T — cos(wt)] , , ^ , 

- 5 -^ coth{Pw/2) 

. [sin{wt) — wt] 

-f 1 5 . 
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Usually we are only interested in the coherence decay, 
and so the imaginary terms in the exponentials above are 
often neglected. As discussed earlier, the coherence decay 
is dominated by the behaviour of J(w) at low frequencies. 
To demonstrate this we characterize the low frequency 
components of the spectral function within the standard 
form J{w) = [II, 0- The behaviour of 

the coherence decay depends critically on the parameter 
s. At any finite temperature, if 0 < s < 2 then the 
coherence |pi 2 (i)i will decay exponentially to zero in the 
long time limit 0 . The rate of this decay is linear in the 
spectral coefficient a, and the decay becomes arbitrarily 
slow as s approaches 2 00- Meanwhile if s > 2 then 
the coherence decays initially, before settling at a finite 
value |pi 2 (oo)| >0 0 . The evolution of coherence for 
s = 1, 2 and 3 is exhibited in figure [2l Consequently, we 
might hope to suppress the decoherence of our subsystem 
by engineering a spectral density which suppresses the 
growth of the spectral density at low frequency; such 
that the coefficient a corresponding to any component 
with s < 2 is small. This idea has been proposed in a 
range of contexts, most recently as a means of sustaining 
coherences within photosynthetic systems 0. 

The analysis above has considered the spectral density 
to be a single well defined function. In practice, most 
experiments are performed over an ensemble of systems. 
Additionally in the photosynthetic context, excited states 
coherences within a single system must be preserved be¬ 
tween the different molecular sites contributing to trans¬ 
port. Each site will be coupled to a different vibrational 


bath. Since any ensemble will exhibit a range of spectral 
densities; we must quantify the effect spectral disorder 
has on the coherence poiit)- 

In the above we neglected the oscillatory terms. For 
later convenience we note that, in the absence of spectral 
disorder, the complex term oscillates as Exp[—i{Ei — 
ar(s)i(;c)i] in the long time limit. We may identify the 
reorganization energy A = Q!r(s)r(;c. We also note that 

r(3) = 2. 



t 


FIG. 2: Coherence decay for spectral functions with s = 1, 
2 and 3. The reorganization energy A is fixed to 0.1, and 
we have chosen the energy scale Wc = 1. We have used the 
high temperature limit coth{l5w/2) 2kBT/w, and chosen 
keT = 10 . 


III. DISORDER IN THE SPECTRAL 
COEFFICIENT 


For the remainder of this paper we consider the special 
case s = 3, which is a common form of spectral function 
arising from the coupling to bulk phonons 0, and we 
choose our energy scale such that Wc = I- Thus the 
mean spectral density J(w) = aw^e~'^. Furthermore, 
we enter the high temperature limit Pwc <C 1. In this 
limit we may expand coth{(dw/2) —>■ 2kBT/w. Recalling 
equations im we define 


ip{t) = - f dwJ{w)f{w,t) 
a Jo 


( 7 ) 

( 8 ) 


Under our approximations 
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The simplest way to introduce disorder to our system is 
to allow the spectral coefficient a to vary, while holding 
the spectral distribution constant. This situation, where 
the reorganization energies of different, spatially sepa¬ 
rated pigments varies was proposed for pigment-protein 
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complexes in Refs. To do this, we average across 

the probability distribution of a, denoted by P{x). 

Pi2{t) = ( 11 ) 

poo 

= / dxP{x)e-^^^^\ ( 12 ) 

Jo 

The coherence evolution will depend on the choice of 
probability distribution P{x). For simplicity, we assume 
that P{x) is Gaussian distributed about the mean value 
a with standard deviation a. 

poo 

(j)(t) = f{a,a) (13) 

Jo 

The prefactor f{a,a) ensures normalization. Equation 
M has an exact solution in terms of the error function. 
This solution is used in figure [3l where we plot the co¬ 
herence oscillation for fixed ct at a range of temperatures. 
To understand the behaviour physically we assume that 
CT <C a and extend the lower limit of integration to — c» 
to obtain 

= Exp[-aipR - cr^((/5? - <Pfl)/2] X 

Exp[-i{a-a'^(pR)(pi]. (15) 


Extending the lower limit of integration requires a^cpR ^ 
a; this ensures that the real part of /(w) induces a coher¬ 
ence decay as expected. Thus the result above can only 
be valid for spectral densities with s > 2, since PR{t) 
must be bounded. Neglecting this small term we obtain 
the coherence decay 


|/Ooi(i)| = Exp[-ap{t) - 


= Exp 


2akBTt^ 

l + P 



1 

( 1 +^ 


(16) 

(17) 


In the absence of disorder the coherence stabilizes at a 
minimum value, |pi 2 (oo)| = Exp{—2akBT). Introduc¬ 
ing disorder causes the coherence to decay to zero on a 
timescale l/cr. In hindsight, the origin of this decay is 
clear. The disorder in the spectral function causes the re¬ 
organization energy to become disordered, at times t > 1 
this introduces a static energetic disorder into the tran¬ 
sition energy, which drives a coherence decay over the 
ensemble. This dephasing term is suppressed for t < 1, 
since the vibrational bath takes a finite time to relax. 
In the case considered, the reorganization energy is also 
Gaussian distributed with standard deviation a a = 2a. 

More striking however, is the oscillating term in Eq. 
[Thl Denoting ujBt = (a — a'^ipR)(pi, and entering the 
long time limit t ^ 1 we find 

OJB = -2{a-2a^kBT) (18) 

= -A + Aa^kBT (19) 

The disorder has introduced a small shift in the oscil¬ 
lation frequency of the coherence, which grows linearly 


as temperature increases. Under our earlier approxima¬ 
tions, this expression fails when 2a^kBT ~ a. 

This shift has an intuitive origin. When we introduce 
disorder to our spectral function, the variation in the real 
coherence decay (driven by pr) is correlated to variation 
in the reorganization energy A. Gonsequently the aver¬ 
age over oscillation frequencies becomes a weighted aver¬ 
age favouring higher frequencies (smaller reorganization 
energies), and the apparent oscillation frequency mea¬ 
sured experimentally differs from the true mean value. 
The shift grows with temperature, since pR oc ksT. 
Reinstating Wc we find that dujB/dT = a\kBlw\, or 
0.7 ((ta/wc)^ cm“^/K. a A labels the standard deviation 
of the reorganization energy. We note that wb is inde¬ 
pendent of time when t ^ 1 because PR{t) is bounded, 
this would not be the case if s < 2. Although this shift 
is small, it could potentially be observed experimentally 
in appropriate systems. Though we have considered a 
specific spectral function here, a similar effect may arise 
whenever a system exhibits significant disorder in the 
density and coupling strength of low frequency modes. 

Indeed such an effect could underlie the increase of 
beating frequency with temperature observed by Pan- 
itchayangkoon et al [T^ . Panitchayangkoon tentatively 
observed a shift of ^5 cm“^ between 77 and 150K, in the 
coherence between excitons 1 and 3 of the EMO complex. 
Accounting for such a shift within our model requires 
that CTA 0.3u>c = O.lwp, where Wp labels the frequency 
of the peak in the spectral function. There remains a 
great deal of uncertainty regarding the spectral function 
of the EMO complex , and to our knowledge there is 
no published work on the statistics of spectral function 
disorder within an ensemble. Most studies have not used 
the simple cubic spectral function considered here. How¬ 
ever, as a rough estimate we note that the smooth, low 
frequency part of these spectral functions typically have 
values for Wp in the range 150-250 cm“^. This implies 
a standard deviation a a 15-25 cm“^. Since typical 
values for A lie in the range 20-60 cm“^ [^, [T3, [Ta. [^. 
this suggests that the ratio cta/A ~ 0.25-1.25. Given the 
limited data available, we conclude that spectral function 
disorder could plausibly generate a shift on the same or¬ 
der of magnitude as that observed experimentally, though 
many other factors would need to be assessed in order to 
make any definite assignment. 


IV. DISORDER IN THE SPECTRAL 
DISTRIBUTION 

The previous section illustrated the key effects of dis¬ 
order on the coherence evolution within a simple model. 
Here we demonstrate that the same behaviour is observed 
when a more realistic continuously varying spectral dis¬ 
order is considered. 

We consider the spectral function J{uj) = J{uj)+aj{uj). 
aj{uj) represents some random disorder at the frequency 
w, and is correlated over a frequency range ^ -C 1. Aver- 
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FIG. 3: (A) Coherence oscillations at a range of temperatures. 
We fix Wc = Wp/S = 50 cm“^, A = 40 cm“'^, cta = 15 
cm“^, and the underlying energy gap Ei = 200 cm“^. For the 
dashed curve, the disorder a a = 0. (B/C) We fit these curves 
to the function cos{u)t)e and extract estimates for the 
frequency shift and dephasing rate. We predict the frequency 
shift using eqn. 19, while the dephasing rate is described by 
a line of best fit. We note that for the parameters considered 
here, the frequency shift is similar to that observed in Ref. 
but the dephasing rate is roughly twice as large. 


aging over disorder realizations implies 

pi2(t) = (20) 

poo 

(j){t) = {Exp[ dwaj{w)f{w)]). (21) 
do 

As previously we take J{w) = aw^e~'^, and for simplicity 
we take <Jj{w) = J{w)a(w)/a. We will assume that the 
probability distribution for cr is symmetric about 0, which 
implies that <j{w) is a random number bounded between 
±a. We may express (f){t) as 

/ OO pOO 

/ . ( 22 ) 

-OO —OO 


P(a:, y) represents the probability distribution of the two 
correlated variables 


X = 


dwe ™(1 — cos{wt))a{w), (23) 

Jo 

poo 

y= dwe~'“ w{sin{wt) — wt)a{w). (24) 

Jo 


At this point we note that both x and y have zero mean. 
Furthermore both are effectively described by a sum over 
0(1/^) uncorrelated random components. Consequently 


so long as ^ is small we may apply the central limit the¬ 
orem and describe P{x, y) via the normal distribution of 
two correlated variables, 


P{x,y) = 


Exp[- 


27rcr3,cry Y^l - 


2^V 


2(1-p^) cr2 a^^ayj 


(25) 


In the expression above, cr^ =< x^ >, =< y^ >, and 

p = OxylioxCJy) where a^y =< xy >. Inserting P{x,y) 
into equation 1221 we find 




(26) 


Note that in equation [22] we extended the limits of the 
integrals over x and y to ±oo, this is only valid if cr^ <S 
a/ksT. In the appendix we show that 


<Xxy{t) = (3 ^ 


poo 

/ dwe~‘^'^{l — cos{wt))^, (27) 

Jo 

pOO 

/ dwe~^'^w^{sin{wt) — wt)^, (28) 

Jo 

poo 

/ dwe~'^'“w X 

Jo 

(1 — cos{wt)){sin[wt) — wt). (29) 


The random variable a[w) has been extracted and re¬ 
placed by its variance /3^. Additionally, all three vari¬ 
ances are proportional to the correlation length Since 
tr is bounded by ±a, jd < a. The exact solution of these 
three integrals is given in the appendix, here we simply 
note that in the long time limit t ^ I, cr^ = 3/3^^/4, 
(Jy{t) = 3f^/3^^/4, and (Jxyit) = —t/3^C/4. Thus at long 
times 

(j){t) Ki Exp[ --— \Exp[-il3‘^^kBTt]. (30) 

O 


We insert equation |30| into equation |20| to describe the 
coherence oscillations at long times. Just like in the pre¬ 
vious section, the spectral disorder drives the coherence 
amplitude to zero in the long time limit, and a tem¬ 
perature dependent shift in the apparent oscillation fre¬ 
quency of the coherence is observed. As before, we may 
obtain the standard deviation of the reorganization en¬ 
ergy, a A = /3v3j/4, S'lid the shift in the coherence fre¬ 
quency, ujB = —A -I- fd^^kBT. Reinserting Wc, we find 
dujB/dP = Aa\kB/‘3wl. 


V. CONCLUSIONS 

In this work we investigated the effect of disorder in the 
vibrational spectral function of an ensemble of two level 
systems. We find that disorder introduces an additional 
coherence decay, since it introduces disorder in the reor¬ 
ganization energy of the excited state. This additional 
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source of decoherence is suppressed at early times, as the 
vibrational bath exhibits a finite relaxation timescale. As 
well as this coherence decay, the apparent oscillation fre¬ 
quency of the coherence is slightly increased. This shift 
occurs because the decoherence of a single realization of 
the system is correlated to its reorganization energy. The 
shift increases with temperature, an effect which may be 
observed experimentally in ultrafast spectroscopies [l9l |. 


y K, ^ e ^^n^{sin{n^t) — n^t)(Tn. 

n 

Note that we have explicitly imposed correlation over a 
frequency scale ^; an now refers to independent uncorre¬ 
lated random numbers. All these random numbers share 
a common variance < tr^ >= /3^. As an example we 
consider the covariance < xy >, this is given by 
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Appendix 

We wish to find the variances and covariance of 

pOO 

x= dwe~'^{l — cos{wt))a{w), 

Jo 

poo 

y= dwe~'^w{sin{wt) — wt)a{w). 

Jo 

a(w) represents a random function correlated over an 
interval Since the correlation length ^ 1 is short, 

we may approximate 

X K, ^ e~"^(l — cos{n^t))an, 


< xy > ~ ^ e — cos{n^t)){sin{n^t) — n^t) 

n 

poo 

= 0^^ / dwe~^'^w{l — cos{wt)){sin{wt) — wt). 

Jo 

Following this procedure 

poo 

< — cos{wt))‘^, 

Jo 

pOO 

< y^ > = / dwe~^'^{sin{wt) — wt0. 

Jo 

Performing the three integrals, we obtain 


< a;" > = 


< 2 / > = 


(„ 768{f - 4) 6 + 3f + 


6 


(^2+4)4 (t2 + l)3 


-b 106t2 -b 14t4 + t6) 
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